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We present results for two-particle transverse momentum correlations, {Ap t ,iAp t ,j}, as a function 
of event centrality for Au+Au collisions at ^/sjvjv = 20, 62, 130, and 200 GeV at the Relativistic 
Heavy Ion Collider. We observe correlations decreasing with centrality that are similar at all four 
incident energies. The correlations multiplied by the multiplicity density increase with incident 
energy and the centrality dependence may show evidence of processes such as thermalization, jet 
production, or the saturation of transverse flow. The square root of the correlations divided by the 
event-wise average transverse momentum per event shows little or no beam energy dependence and 
generally agrees with previous measurements at the Super Proton Synchrotron. 

PACS numbers: 25.75.Gz 



The study of event-by-event fluctuations in global in particle production, may provide evidence for the pro- 
quantities, which arc intimately related to correlations duction of the quark gluon plasma (QGP) in relativistic 
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heavy ion collisions [JllllSll@,EBIl[ll[lJ[ll 

EJ E3- Various theoretical work predicts that the 
production of a QGP phase in relativistic heavy ion col- 
lisions could produce significant dynamic event-by-event 
fluctuations in apparent temperature, mean transverse 
momentum, multiplicity, and conserved quantities such 
as net charge. Several recent experimental studies at the 
SPS 0E113 and at RHIC Ea E3, Hl|, H3, H3, H3| have 
focused on the study of fluctuations and correlations in 
relativistic heavy ion collisions. One possible signal of 
the QGP would be a non-monotonic change in p t corre- 
lations as function of ccntrality and/or as the incident 
energy is raised §]. 

Here we report an experimental study of the incident 
energy dependence of p t correlations using Au+Au colli- 
sions ranging in center of mass energy from the highest 
SPS energy to the highest RHIC energy using the STAR 
detector at RHIC. 

Fluctuations involve a purely statistical component 
arising from the stochastic nature of particle production 
and detection processes, as well as a dynamic component 
determined by correlations arising in various particle pro- 
duction processes. In this paper we first unambiguously 
demonstrate the existence of a finite dynamical compo- 
nent at all four incident energies by comparing the dis- 
tribution of measured event-wise average transverse mo- 
mentum per event, (pt), with the same quantity from 
mixed events. We then analyze these dynamical fluctua- 
tions using the two particle transverse momentum corre- 
lations defined as covariance 
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and A event is the number of events, pt^ is the transverse 
momentum of the i th track in each event, A/- is the num- 
ber of tracks in the k th event. The overall event average 
transverse momentum ((pt)) is given by 
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where (p*) fc is the average transverse momentum per 
event given by 



FIG. 1: (Color online) Histograms of the average transverse 
momentum per event for Au+Au at ^snn = 20, 62, 130, 
and 200 GeV for the 5% most central collisions at each en- 
ergy. Both data and mixed events are shown for each incident 
energy. The lines represent gamma distributions. 
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(Apt t iAp t j) is independent, to first order, of detection 
efficiencies because both the numerator Ck and the de- 
nominator Afc(Afc — 1) are proportional to the the square 
of the particle detection efficiency. Therefore the ef- 
ficiency cancels. By construction (Apt. i Apt. j) is zero 



within statistics for properly mixed events because all 
correlations are removed. Note that we use mixed events 
only in Figure 1. 

The data used in this analysis were measured using the 
Solenoidal Tracker at RHIC (STAR) detector to study 
Au+Au collisions at ^/s]^ = 20, 62, 130, and 200 GeV 
PH . The main detector was the Time Projection Cham- 
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ber (TPC) located in a solenoidal magnetic field. The 
magnetic field was 0.25 T for the 20 and 130 GeV data 
and 0.5 T for the 62 and 200 GeV data. Tracks from the 
TPC with 0.15 GeV/c < p t < 2.0 GeV/c with |r?| < 1.0 
were used in the analysis. All tracks were required to 
have originated within 1 cm of the measured event ver- 
tex. Events were selected according to their distance of 
the event vertex from the center of STAR. Events were 
accepted within 1 cm of the center of STAR in the plane 
perpendicular to the beam direction. For the 20 and 
130 GeV data sets, events were accepted with vertices 
within 75 cm of the center of STAR in the beam direc- 
tion, while for the 62 and 200 GeV data sets, events were 
accepted within 25 cm of the center. 

Data shown for 62, 130 and 200 GeV are from min- 
imum bias triggers. Minimum bias triggers were de- 
fined by the coincidence of two Zero Degree Calorimeters 
(ZDCs) [26j located i 18 m from the center of the in- 
teraction region. For 20 GeV a combination of minimum 
bias and central triggers was used. Centrality bins were 
determined using the multiplicity of all charged particles 
measured in the TPC with \rj\ < 0.5. The centrality bins 
were calculated as fractions of this multiplicity distribu- 
tion starting with the highest multiplicities. The ranges 
used were 0-5% (most central), 5-10%, 10-20%, 20-30%, 
30-40%, 40-50%, 50-60%, 60-70%, and 70-80% (most pe- 
ripheral) . Each centrality was associated with a number 



A, 



part 



using a Glauber Monte 



of participating nuclcons, 
Carlo calculation [2tT| . 

We treated the variation of ((pt)) within a given cen- 
trality bin using the following procedure. We calculated 
((pt)} as a function of N c h, the multiplicity used to define 
the centrality bin. We fitted this dependence and used 
the fit in Eqs. 1-4 on an event-by-event basis as a func- 
tion of N c h ■ This method removes the dependence of the 
experimental results on the size of the centrality bin and 
slightly reduces (Ap t ^Ap t j) by removing correlations in- 
duced by the changing of ((pt)) within the experimental 
centrality bins. The results presented in the paper are 
obtained using this fitting procedure. 

Fig. shows histograms of (p t ) for the 5% most cen- 
tral Au+Au collisions at 20, 62, 130, and 200 GeV. His- 
tograms for (p t ) are also shown for mixed events. The 
histograms for the data are wider than the histograms for 
mixed events indicating that we observe non-statistical 
fluctuations at all four incident energies. Similar results 
are obtained for all centralities. The overall normaliza- 
tion reflects the number events taken at each energy. The 
values of pt included in these histograms are not corrected 
for experimental momentum resolution, acceptance or ef- 
ficiency. 

The mixed events at each energy were created by ran- 
domly selecting one track from an event chosen from 
measured events in the same centrality and event ver- 
tex bin. Ten centrality bins and either five or ten bins 
(depending on the available number of events at each en- 
ergy) in the event vertex position in the beam direction 
were used to create mixed events with the same multi- 
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FIG. 2: (Color online) (Apt,iApt,j) as a function of central- 
ity and incident energy for Au+Au collisions compared with 
HIJING results. 



plicity distribution as the real events. Note that we do 
not use mixed events for the quantitative analysis based 
on (Ap t ,iAp t ,j). 

The lines in Fig. ^ represent gamma distributions for 
both the data and mixed events. The parameters for the 
gamma distributions are shown in Table 1. According 
to Ref. |28|. without pt cuts the parameter a divided by 
the average multiplicity in the centrality bin, (A), should 
be approximately two and the parameter [3 multiplied by 
(N) should reflect the temperature parameter of the pt 
distributions. We find that a/(N) varies from 2.27 to 
1.93 and (3(N) varies from 0.230 to 0.299 GeV/c as the 
energy goes from 20 to 200 GeV. 



TABLE I: Parameters for the gamma distributions shown in 
Fig. 1. The gamma distribution is given by the form f(x) = 

- r^ffa 13 where a — ^ and (5 = — in GeV/c. p, is the 
mean in GeV/c and a is the standard deviation in GeV/c. 
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To characterize the transverse momentum correlations, 
we use the the quantity (Ap ti iAp t j} , defined in Eq. 1. 
Fig.Oshows (Apt^Aptj) for Au+Au collisions at ^snn 
= 20, 62, 130, and 200 GeV as a function of centrality. 
One observes that (Ap t i Ap t j) decreases with central- 
ity at all four energies as expected due to a progressive 
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FIG. 3: (Color online) (dN/dri){Apt.iAptj) as a function of 
centrality and incident energy for Au+Au collisions compared 
with HIJING results. 
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FIG. 4: (Color online) (Apt,i Aptj) / {{pt)) as a function of 
centrality and incident energy for Au+Au collisions compared 
with HIJING results for corresponding systems. The inset 
shows the excitation function for the most central bin. 



dilution of the correlations resulting from the increased 
number of participants if the correlations are dominated 
by pairs of particles that originate from the same nucleon- 
nucleon collision. The correlations measured at 62, 130 
and 200 GeV are similar while the correlations for 20 GeV 
are smaller than those observed at the higher energies. 

To explore the issue of the relative importance of short 
range correlations such as Coulomb interactions and Han- 
bury Brown-Twiss (HBT) effects, we extracted the cor- 
relations excluding pairs with invariant relative momen- 
tum, qi nv , less than 0.1 GeV/c, assuming that all parti- 



cles were pions. We observed that 10% of the measured 
correlations at 62, 130, and 200 GeV and 20% of mea- 
sured correlations at 20 GeV could be attributed to these 
short range correlations. These estimates agree with 
those extracted for 17 GeV Pb+Pb 0] using a somewhat 
different method. We also estimated the contribution of 
resonances and other charge-ordering effects by studying 
the reduction in the correlations for same charge (nega- 
tive) particles compared with correlations for all charged 
particles. This study indicated that the reduction in 
(Apt,iAp t ,j) is 40% at 20 GeV, 20% at 62 and 130 GeV, 
and 15% at 200 GeV. We do not correct (Ap t>i Ap t j) for 
short range correlations or resonance contributions. 

The errors shown in all figures are statistical unless 
otherwise noted. We estimate the systematic relative 
errors for (AptjAptj) using studies of the effects of 
Pt-dependent efficiencies (1.2%) and sensitivity to track 
merging and splitting (1.4%). These values give an over- 
all systematic relative error of 2%. The measured cor- 
relations were lowered approximately 3% using the fit- 
ting method rather than the binning method. The re- 
ported values are sensitive to the y± cuts for kinematic 
and physics reasons. Using HIJING [2!j we observe a 6% 
increase in correlations when the lower p t cut is removed. 
Raising the upper p t cut increases the correlations. We 
used 0.15 GeV/c < p t < 2.0 GcV/c for all the results 
reported in this paper. The upper p t cut was chosen to 
be consistent with previous work [l3.l24|. 

Also shown in Fig. [5] are HIJING calculations for 
Au+Au collisions at y/s^N = 20, 62, 130, and 200 GeV 
p9j . We used HIJING version 1.36 with the default op- 
tions, which includes jet quenching. The HIJING results 
were obtained by selecting particles with 0.15 GeV/c 
< Pt < 2.0 GeV/c with \r)\ < 1.0 without further effi- 
ciency corrections. HIJING reproduces correlations in 
p+p and a+a collisions at Intersecting Storage Ring 
(ISR) energies [3(J, p+p collisions at RHIC energies, and 
p+p collisions at CERN p+p Collider (SppS) energies 
|31|. We use HIJING to provide a reference incorporat- 
ing a superposition of nucleon-nucleon interactions. Any 
differences between HIJING and the experimental results 
might signal phenomena unique to nucleus-nucleus colli- 
sions. The HIJING calculations exhibit little incident 
energy dependence and decrease with increasing central- 
ity. The values for (Ap t ,iAp t ,j) predicted by HIJING are 
always smaller than the data. 

To address the observed dilution of the correla- 
tions with centrality, and to check the hypothesis 
that the correlations scale as inverse multiplicity, we 
multiply ( Ap t , i Ap t j ) by the charged particle pseudora- 
pidity density at a given centrality, dN/drj. We use 
fully corrected values for dN/drj from published work 
III HI IP- The quantity ^(Ap M Ap M ) then is in- 
sensitive to efficiency and is similar to the (efficiency 
corrected) quantity Aa pt that STAR has reported 
previously. 

In Fig. we show the quantity (AptjAptj) for 
Au+Au collisions at 20, 62, 130, and 200 GeV as a func- 



6 



tion of centrality. In this figure the errors include the 
quoted errors in dN/drj. This quantity increases with 
incident energy at all centralities. At each energy this 
measure of the correlations increases quickly as the colli- 
sions become more central and then saturates in central 
collisions. The behavior of this quantity is similar to that 
of the quantity Aa pt previously studied by STAR |19| . 
This saturation might indicate effects such as the onset 
of thermalization [13, the onset of jet quenching [l4|, the 
saturation of transverse flow [35| in central collisions, or 
other processes. 

In Fig. |31 the results of HIJING calculations for 
£^-{AptjAptj) are also shown. In contrast to the ex- 
perimental results, the HIJING results show little de- 
pendence on centrality. 

To account for possible changes of (Ap t ^Ap t j) due to 
possible changes in ((pt)) with incident energy and/or 
centrality of the collision, we also study the square root 
of the measured correlations scaled by ((pt))- The re- 
sulting quantity yj (Ap t ^Ap t .j) / {(pt)) is shown in Fig.0] 
for Au+Au collisions at 20, 62, 130, and 200 GeV. Sim- 
ilar results from Pb+Pb collisions at 17 GeV ^(| are 
also shown in Fig. These values are consistent with 
our measured results for Au+Au at 20 GeV. We observe 
little or no dependence on the incident energy for this 
quantity. The inset in Fig. 0] demonstrates the incident 
energy dependence of yj (AptjAptjj/ ((pt)) for the - 
5% most central bin where the Pb+Pb results are from 
Ref. [3. 
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